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Abstract 

We consider the problem of option pricing and hedging when stock returns 
are correlated in time. Within a quadratic-risk minimisation scheme, we 
obtain a general formula, valid for weakly correlated non-Gaussian processes. 
We show that for Gaussian price increments, the correlations are irrelevant, 
and the Black-Scholes formula holds with the volatility of the price increments 
on the scale of the re-hedging. For non-Gaussian processes, further non trivial 
corrections to the 'smile' are brought about by the correlations, even when the 
hedge is the Black-Scholes A-hedge. We introduce a compact notation which 
eases the computations and could be of use to deal with more complicated 
models. 



1 Introduction 

The assumptions underlying the Black-Scholes model [|l], § - that the dynamics of 
financial assets can be modelled by a continuous time Gaussian process - are very 
far from reality. It is well known, in particular, that the empirical distribution of 
returns exhibits 'fat-tails'. The presence of these extreme events not only induces 
a 'smile' in the implied volatility of the options, but also destroys the essence of 
the Black-Scholes pricing procedure: the possibility of constructing a riskless hedge 
which replicates perfectly the option pay-ofF P, |^. A more subtle further effect, 
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completely discarded by the standard Black-Sclioles model, is the presence of small, 
but significant, temporal anticorrelations in the series of stock returns. Correspond- 
ingly, the variance of the A^-day distribution of returns is smaller than times the 
variance of the daily distribution (see e.g. 0). On other financial assets, such as 
stock indices, one observes the opposite effect of weak positive correlations. The or- 
der of magnitude of these correlations or anticorrelations is 10% from one day to the 
next. These correlations lead, in principle, to the possibility of statistical arbitrage. 
However, one should also take into account transaction costs which significantly re- 
duce the possibility of using these correlation in practice [^]. In other words, the 
presence of transaction costs allow the existence of non zero correlations, which can 
be significant and therefore must be dealt with consistently in the context of option 
pricing. 

What is the infiuence of these correlations on the price of an option and on 
the corresponding hedging strategy? The aim of this paper is to provide a general 
answer to this question within a quadratic risk minimization sheme, in the case 
where these correlations are small [| (which is, in most cases of interest, amply 
justified), with only very few assumptions on the actual distribution of returns. 
We calculate the first order corrections both to the price and optimal hedge. In 
the Gaussian case, we find that the effect of correlations can be compensated by 
a change in the hedging strategy and therefore options should be priced using the 
standard uncorrelated Black-Scholes model 0. The correct volatility to be used is 
the one measured on the time scale of the rehedging, and not the one corresponding 
to the terminal distribution - i.e. corresponding to the maturity of the option. 
In other words, the 'risk-neutral' measure that one should use to price derivatives 
is simply an uncorrelated Gaussian measure with the same elementary variance. 
These conclusions are not valid if one considers a process for which the underlying 
distribution is non-Gaussian. One then finds contributions to the option prices 
which depend in general on the strength of the correlations, and on the higher 
moments of the distribution (in particular on the kurtosis). The corrections to 
the option price also depend on the past history of the price-increments of the 
underlying. The detailed analysis of a specific example shows, in particular, that 
anticorrelations change the effective kurtosis to be used when computing 'smile' 
corrections to the options prices. In general, the correct kurtosis to be used is 
neither the local nor the global one, but some effective kurtosis which depends on 
the strength of the correlations. 



^The general correlated Gaussian case was considered in and leads to rather complex for- 
mulae for which no simple interpretation was given. 

^This can be easily seen within a Black-Scholes framework: see Appendix. 
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2 General setting 



Our starting point is the global wealth balance approach discussed in P, ^. We 
set, for simplicity, the interest rate to zero, and consider a discrete time model with 
time interval r, within which no rehedging takes place. The global wealth balance 
at time T = Nt for the writer of the option reads 

W = C + H -Q, 

where C is the option premium, is the pay-off of the option, and H is the result 
of the hedging strategy 

7V-1 

i/=^0„A„+i. (1) 

n=0 

In equation (0), A„+i = Xn+i denotes the change of the price x of the underlying 
asset between time n and n+ 1, and 0„ is the amount of underlying in the portfolio 
at time n 0. 

In incomplete markets, the price of an option is not unique since a risk-premium 
should be added. However, a convenient framework is the risk minimisation scheme 
proposed in |]^, ^, |^, that fixes both the price and strategy C, 0„ so as to minimize 
the total risk TZ associated to option writing, which we simply define as 

7^2 = (iy2^ , 

where (...) is an average over the objective probability at time 0. Minimizing with 
respect to C fixes the price of the call option to be 

C={Q-H). (2) 

If, on the other hand, we minimize TZ^ with respect to 0„, we obtain an equation 
implicitly determining the optimal strategy 0*, equation that reads, for < < 

(A„+i(0-if))„ = C(A„+i)„ (3) 

The notation (...)„ in the above equation means that the average (again over the 
objective probability) is performed at time n - i.e. with all the information available 
at this instant of time. This notation turns out to be very powerful, and could be 
used to handle very general situations. We now analyze the above general and 
compact equations (0) and @ in several simple but important limiting cases. 

■^In the following, we assume for simplicity that the interest rate is zero. It is easy to generalize 
the formulae obtained here when this rate is non-zero: see Istl for details. 
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3 The case of zero conditional drift 

Let us define the conditional drift /i„ and tlie conditional variance D„ at time n as 

l^n = {^n+l)n 
Dn = ((A„+i-/i„)'>„- 

Note that these quantities are the expected drift and variance for the next time 
interval, and depend in general on the past (realized) history of price changes. 

We start by assuming that the drift /i„ is identically zero (this is to say, the price 
process is a 'martingale'). Then, since (pn is known at time n, one has that 

(0„A„+i) = (0„(A„+i)„) = 0, 

and therefore that (H) = 0. This implies that the option price is nothing but the 
expectation value of the pay-off 

C = (9) . 

To determine the optimal hedge, we note that, in general, one has 

(An+iH)^ = ^ (0^A„+iA„+i)„ (4) 

m=0 

n-1 N-1 

= (Pn (^n+l)„ + X] ^"lArn+lfl'n + ^ (0m/imAn+l)„ . 
m=0 m=n+l 

In the case of vanishing drift /i„ = we therefore have that 

(A„+iiy)„ = 

so that the optimal hedging strategy can be computed, from equation (^, as 0, |[ 

K = ^ (An+10)n • 



4 Perturbative expansion for small drift 

Let us now turn to the case where the conditional bias /i„ is non-zero but small. 
Using equation (^), it is easy to see that the basic equations (|]) and (^, which 
determine the optimal price and the hedging strategy, can be solved perturbatively 
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in Hn- For our purposes, we will only need the first correction in the conditional 
drift, which we now derive. First note that equation (Q) can be rewritten as 

N-l 
n=0 

In particular, to first order in one finds that 

c=(e)-$:(^^0). (5) 

n=0 \ " / 

One can also compute the first correction to the hedging strategy, by combining 
equation @ with the zero-th order results obtained in the last section. One obtains 
then 

^S-l^(A.e>-^|'(^A.e)-i-,.(e). (6) 

In the above equation, we have only shown the optimal hedge for the first time 
period 0g. The explicit expression for the optimal hedge for the subsequent time 
periods is also of a similar form, but it is quite intricate and not very illuminating. 
More importantly, it is not necessary to determine the correct hedging strategy. 
This important point will be discussed, in a slightly more general context, in section 

0- 

5 A general model for weakly correlated processes 

Let us introduce a rather general model for correlated price increments. We first 
define an auxiliary set of uncorrelated i.i.d. random variables {rjn}, distributed 
according to an arbitrary probability distribution]^ P{r]). The joint distribution is 
then given by 

P{{Vn})=llP{Vn). 
n 

We will assume that P{ri) is symmetrical^ - i.e. that 

j 7f''+^P{7])d7] = 0, (7) 

^The following results can be generalized to the case where the distribution P{ri) explicitly 
depends on n. 

'"'In the sequel, we will actually only need that equation (1^) holds for g = 0,1. 
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and we will denote the second moment of P by 

D = j TfP{r])dT]. 

We now construct the set of correlated price increments {A„} by writing 

A„ = ^ MnmVm + (8) 
m 

with 

The coefficients c„ are assumed to satisfy 

Co = c-k = Ck- 

In the sequel, we assume that both the c^'s and fi are small, and will work only to 
first order in both fj, and c. 

Let us show the significance of equation (H), by first noting that, since {rjn) = 0, 
we have that 

(A„) = fi. 

Therefore /x is nothing but the average unconditional drift of the price process. 
Moreover, to first order in c, the following holds 

((A„ - fi){Am - Ai)) = MniMmj {ViVj) = D5n,m + C„-m- 

The coefficients D and c denote then, respectively, the unconditional variance and 
correlation of the price increments. Finally, one has that, independently of c„-m, 
and for all p, 

((A„-/xr) = (^o- 

We note that the variance of the price differences over a time scale n, given by 

/ / m+n \ 2\ 

\ \fc=m+l / / 

is related to the correlation coefficients as follows 

n-l 

a'^{n) = nD + ^2{n-k)ck 

k=l 

Cn = ^{a\n + l)-2a\n) + a^{n-l)). (9) 
6 



In the case where c„ = 0, one finds a'^{n) = nD, whereas, in the case of very 
short range anticorrelations (c„ = —c6n,i), one has (j^(n) = n{D — 2c) + 2c. In this 
case, the global volatility a'^{n)/n D — 2c is lower than the short-range volatility 
a\l) = D. 

The probability distribution of the price increments can be written, to first order 
in /i, c, as 



dAn 2D 



m / 

1 ^d\nP 



(10) 



c A 



It is now easy to write the marginal distribution at time n = 0, given all previous 
returns A„ with n < 



E 

n>0 



ainP 



2D ^ 



51nP 



max(n,m)>0 



9A, 



-^n—m'—^m 



A. 



n>0 



The equation above follows immediately from the general expression for the prob- 
ability distribution, up to a possible normalization factor. To show that (|1T|) is 
actually correctly normalized, we have only to note the obvious fact that 



dhiP 



Moreover, using that 



ciAP(A)A 
rfAP(A)A' 



(91nP 

dA 
.d\nP 
dA 



-1 




it is then easy to show that the conditional drift and variance /i„ and Dn are given, 
to first order in fi and c, by 



/i 



D. 



l91nP 
2 SAT 



6 Derivative pricing with small correlations 
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6.1 A general formula 

Let us now apply the general perturbative formula @ for the price of the contract 
to the specific model of correlations described in the previous section . One finds 



(0) 



N , TV 

D 



n=l n=l m.<n ^ ^ ^ 



Now, the above averaging (■ ■ ■ ) is performed over the 'correlated' historical price in- 
crement distribution, equation (0). It is convenient to re-express all of the following 
formulae in terms of an auxihary probability distribution ]^^^qP(A„), which de- 
scribes an unbiased and uncorrelated process, and is obtained by setting Cfc = yU = 0. 
We will denote the corresponding averages by (...)°, where the superscript indicates 
that the conditional drift is now set to zero (in the case of a risk-free Gaussian model, 
this auxiliary probability distribution is called the equivalent martingale measure in 



the financial mathematics literature). Equation (11) can then be reinterpreted as 



relating the two averaging procedures (■ ■ ■ ) and (■ ■ ■ )° as follows 
, , , ,0 Y-/ 51nP\° 1 >^ / d\nP 

^■■■^ = ^---^ -^^y-^l -2D ^ v"W"-'"^^ 

n>0 ^ " ' max(n,m)>0 ^ 

After a some algebra we then find that 

N , N 



C = (6)° + /i 5^ (F (AO 0)° + ^ E E A„) A^e)° , (12) 

where 



2D 

n=l n=l m<n 



d\nP 1 
F(A) = _ ^A. 

^ ' OA D 

The above equation is the central result of this paper, that we now comment in 
various limits. 

Let us first consider the purely Gaussian case, where InP(A) = —A'^/2D — 
l/21n(27rD). In this situation, it is easy to check that F{A) = 0. Therefore, all 
corrections brought about by the drift and correlations strictly vanish, and the price 
can be calculated by assuming that the process is correlation-free. As shown in the 
Appendix, this is indeed expected in a continuous time framework. In other words, 
the price of the option is given by the Black-Scholes price, with a volatility given 
by the small scale volatility D (that corresponds to the hedging time scale), and 
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not with the volatihty corresponding to the terminal distribution, a'^{N)/N. As 
we have already shown, when price increments are anticorrelated, D is larger than 
a'^{N)/N, and the price of the contract is higher than the objective average pay-off. 
This is due to the fact that the hedging strategy is on average losing money because 
of the anticorrelations: a rise of the price of the underlying leads to an increase of 
the hedge, which is followed (on average) by a drop of the price. 

For an uncorrelated non-Gaussian process, the effect of a non-zero drift fi does 
not vanish, a situation already discussed in details in ||^. The same is true with 
correlations: the correction no longer disappears and one cannot price the option by 
considering a process with the same statistics of price increments but no correlations. 
Note that the correction term involves all with m < 0. This means, as expected 
a priori, that, in the presence of correlations, the price of the option does not only 
depend on the current price of the underlying, but also on all past price increments. 



6.2 The case of short range correlations 

Let us illustrate the general formula in the case of very short range anticorrelations 
Cn = —c6n,i- We will furthermore consider the simple case of path-independent 
pay-offs 0, which are functions only on the terminal price Xn oi the underlying. 
The calculation of the price will therefore involve the following quantities 

h = {F (A„) |x;v)° h = (A„F (A^) |x,v)°, 

where the notation (...|a;7v)'^ means that we are conditioning on a given value of the 
terminal price x^- The computation of the two quantities /i,/2 is easily done in 
momentum space, with the aid the following formulae 

j dA AP{A)e''^ = iDz - ^D^ (3 + k) + o (z^) 

j dA F (A) P{A) e''^ = '-kDz^ + o{z''), 

where k is the kurtosis of the distribution of elementary increments. The result 
reads, for large N, 

K, 4:K 

where pi and p2 are, respectively, the skewness and kurtosis polynomials 

Pi = 7 (x^ - 3x) P2 = 777 (x^ - Qx"^ + 3) 

o 24 

and where we have consider the natural scaling aP' = {xn — a^o)^ /DN : x is the 
moneyness of the option counted in standard deviations. The general equation (|12D 
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then takes the form C = (©)", with a modified pay-off function 

Q[xn) = B{Xn) 1 ( - ) Pi + -7^P2 



y/DN \ 2L)iVy" DN' 

where the exphcit dependence on the last past increment Aq appears. Note that the 
correction terms vanish for at the money options [xn = xo). Let us concentrate on 
the last term in the above equation. It represents an increase in the global kurtosis 
n/N oi the 'risk neutral' distribution that must be used to price the option: 

This is means that the 'implied' kurtosis is increased, in percent terms, as much as 
the global variance is decreased from the local one by correlation effects. There- 
fore the volatility 'smile', which is proportional to the kurtosis, is enhanced by the 
presence of anti-correlations, and reduced by the presence of correlations. Note 
that the above model neglects all volatility fluctuations, which lead to a terminal 
kurtosis decaying much more slowly than n/N (see P]). The above calculations 
could in principle be extended to deal with this effect, which is crucial for flnancial 
applications. 



6.3 An alternative model of correlations 



Let us comment more on the above result. To this end, let us briefly discuss an 
other possible model of correlated price increments. One could consider, instead of 



the probability distribution (TU), the following distribution for price increments (we 
set for simplicity /i = 0) 



(14) 



The above measure induces statistics similar to (|ToD, and in particular one still has 
that (A„Am) = D5n,m + c„_m, and that (A^) = (r/^). On the other hand, it is 
not hard to show, following an analysis similar to the one in section ^, that, to 
flrst order in the correlations, the price of an option is unchanged, regardless of the 
basic day-to-day probability distribution P{i]). How can one decide which model 
more accurately describes a speciflc time-series of price increments? To answer this 
question, we flrst note that (|1^) and (|1^ coincide when P (rj) is Gaussian. It must 
then be that the deviations involve the higher moments of P. Let us denote with Aj 
the i-th cumulant of P (in particular X2 = D and A4 = D^n, with n the kurtosis), 
and with Aj (n) the i-th cumulant of the return Ai + ■ ■ ■ + A„ of the stock in n 
time periods. If the price increments where independent, then Aj (n) = nA,. In the 
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presence of correlations, this equation is not valid any more, but, if the correlations 
are short-range, we can define global cumulants as 



Xi = -\i in) . 



{n oo) 



n 



To be definite, let us consider the case already discussed of short-range anticorrela- 
tions Ci = — c. We have already shown that A2 {n) = n{D — 2c) + 2c and therefore 
that 



It is then not hard to show that, if one considers the model ([10|), then the global 
forth moment is given by 



In both cases, the global historical kurtosis is reduced by anti-correlation effects, 
but the amount varies in the two models. Recalling that 2c/ D ~ 0.10, we see that 
the two models differ by 10% in the measurement of the global kurtosis X^/D^k. 

6.4 The case of Delta-hedging 

Le us conclude this section by discussing other possible hedging schemes. Instead of 
choosing the optimal hedge that minimizes the variance, one can follow the (sub- 
optimal) Black-Scholes A-hedge. As emphasized in this leads in general to a 
negligible risk increase, and has the advantage of removing the effect of the drift n 
on the price, and therefore reduces the risk associated to change of long time trends, 
encoded in /x. One might therefore expect that A-hedging also allows one to get 
rid of the correlations. However, the following calculation shows that this is not the 
case. 

Up to first order in c, it is sufficient to consider that the A-hedge is given by 





On the other hand, if one uses (|14|) one finds that 




d 



(6) 







dlnP 
dA 



(A„+i)e) 



'BS,n — 



dx, 



n 



n 
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The price of the contract is then given by 

N-l 
n=0 

- (e)° + ^f:j:c_(F(AO^e)°. 

n=l m<n 

Note that, in this case, the drift fj, indeed disappears to first order even if F(A) is 
not zero; the correlations, on the other hand, are still present within a A-hedging 
scheme. 



7 Other optimization schemes 

Let us now consider a slightly different problem, where the price of the option is fixed 
by the market to a certain value Cm- The option trader knows that price increments 
have a non-zero conditional drift and wishes to buy and hedge the option in order 
to optimize a certain risk adjusted return. The wealth balance associated to selling 
the option now reads 

w ^Cm + h -e, 

so that the expected return is equal to 

M = Cm + {H-Q)=Cm-C 

and the expected risk is 

We will assume that the trader wants to maximize a certain function G of M and 
71, for example the Sharpe ratio Gi — M/TZ, or a certain risk corrected return 
G2 = XM- n. 

Consider a small variation Scf)^ of the strategy 4>n- The corresponding changes 
in M and TZ read 

— n— = {K+l{H -Q))n+C|ln■ 

We may now extremize G with respect to the hedging strategy and obtain the 
following equation 

8J^dG_ STZdG _ 
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or 

(A„+i(e-i7))„ = /i„c + /i„6; g = n^. (15) 

In the two examples considered before, one has that 

= -- g. = -xn. 

Equation ( ]T5| ) is a generahzation of equation and is again well suited for a 
perturbation in fi. We will work as always to first order in fi, and we will call 
the solution to (^) with the function Q set to zero - that is, the solution to the 
problem solved in section It is then easy to see, using equation (^), that the 
optimal strategy is given by 

T^Sh. (16) 



D. 



n 



We conclude this section by showing that the more general formalism developed 
above actually has some interesting consequences for the original problem of option 
pricing and hedging which was considered in the rest of the paper. Let us in par- 
ticular consider, as we have done in the previous sections, a specific contract sold 
at time at the optimal price Cq and with optimized hedging strategy 0* g (we will 
explicitly insert the time-index for the date of sale of the contract in the last part 
of this section). More generally, we may sell the same contract at some other time 
P at a price Cp with a hedging strategy 0* p (with n> P). It is a natural question 
to ask what is the relation of the various hedging strategies 0* p for fixed time n, as 
we vary the sale time of the contract P. In particular, one may ask if the optimal 
hedging strategy 0* p at time n should be equal to the optimal hedge 0* „ which 
one is to adopt if one is selling the contract at time n. This question can be easily 
answered in the framework of the present section, by noting that, at time n, the 
option traders wealth balance is given by the original price of the contract Cp plus 
the returns of the hedging strategy up to time n 

Wn = Cp + 0p,pAp+i + ■ ■ ■ + 0^_i,pA„. 

From the point of view of the trader, the above quantity acts as an effective market 
price Cm, and the optimization problem after time n is of the type described at 
the beginning of this section. Recalling that the trader is trying to minimize (W"^), 
we deduce that the function G is given by + 7^^ and therefore that Q = M = 
Cm — C = Wn — Cn- We can then use equation (|T6D to deduce that 

Cp = 0n,n - (>Vn - Cn) ■ 
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The above equation has a clear interpretation. In a risk-free model, the quantity 
(W„ — Cn) vanishes, since we can, by assumption, perfectly reproduce the price- 
process of an option by hedging correctly. Therefore the correct hedge is independent 
of the time of sale of the contract. On the other hand, in a model with non-zero 
risk, one needs to correct the hedging strategy whenever the past hedge has not 
completely compensated the price change of the option. In particular, if the hedging 
strategy has made too much money and W„ > C^has to decrease the amount of 
stock in the portfolio whenever the conditional drift yU„ is positive, and increase it 
otherwise. A similar result was obtained in section 3.4; it was however noted 
there that this history dependent strategy, although reducing the variance, actually 
increases the probabilty of large losses. 

8 Summary and conclusion 

We have considered the problem of option pricing and hedging when stock returns 
are correlated in time. Using a variance minimization framework, we have obtained 
a general formula, valid for weakly correlated non-Gaussian processes. We have 
shown that in the limit of Gaussian price increments, the correlations are irrele- 
vant, in the sense that the Black-Scholes formula holds with volatility that of the 
price increments on the time-scale of the re-hedging. For non-Gaussian processes, 
however, further non trivial corrections to the 'smile' effect are brought about by 
the correlations, even when the hedge is the Black-Scholes A-hedge. The kurtosis 
to be used in option-pricing is neither the local one, nor the global one, but some 
effective kurtosis which depends on the strength of the correlations. The above for- 
malism can be extended to treat the case where the distribution of price increments 
depends explicitly on time, as would be needed to treat the case of a time dependent 
volatility. Finally, we have given formula for the optimal hedge when the expected 
return associated to options is non zero. 

Appendix: Correlations within the Black-Scholes 
framework 

The presence of correlations in the continuous time limit can be modelled as a 
history dependent drift term in the stochastic differential equation for the price (or 
log-price process). More precisely, the dynamical evolution of the price reads: 

dX{t) = fx{t)dt + aodW{t) (17) 

where the drift fi{t) depends on the whole past history and dW is a Brownian noise 
of unit variance per unit time. A simple model of correlations is to assume that fi{t) 
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is constructed from the past values of dW{t') as an exponential moving average: 



When e > 0, this corresponds to correlations in price increments extending over a 
time scale F"^. Conversely, e < corresponds to anticorrelations. It is easy to see 
that in this model, the time dependent square volatility cr^(t) is given by: 



that interpolates between the short term volatility ctq and the long term volatility 
Wo + e|- 

In the Black-Scholes formalism, however, the drift is totally absent from the 
price and hedge of an option. This is intimately related to the fact that the Ito 
formula does not involve the drift 0. Hence, it is immediate that in this case, 
correlations are irrelevant and only the short term volatility ctq is needed to price the 
option. This is what we find within our general formalism in the case of Gaussian 
increments. 

Note that the difference between the short term volatility and long term volatility 
is particularly striking in the case of an Ornstein-Uhlenbeck (mean reverting) process 
for which fi{t) = —KX{t). In this case, the long term volatility tends to zero: X{t) 
has bounded fluctuations, even when t ^ oo. Correspondingly, the average pay-off 
over the objective probability distribution is finite for t oo, whereas the option 
price, given by the standard Black-Scholes formula with a volatility (Tq, tends to 
infinity. The difference between the two comes from the hedging strategy that loses, 
because of the anticorrelations, an infinite amount of money in the limit t ^ oo. 
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